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Abstract 

Bianchi type III string cosmological models with bulk viscous fluid for 
massive string are investigated. To get the determinate model of the uni- 
verse, we have assumed that the coefficient of bulk viscosity £ is inversely 
proportional to the expansion 9 in the model and expansion 8 in the model 
is proportional to the shear a. This leads to B = £C n , where £ and n are 
constants. The behaviour of the model in presence and absence of bulk 
viscosity, is discussed. The physical implications of the models are also 
discussed in detail. 
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1 Introduction 



In recent years, there has been considerable interest in string cosmology. Cosmic 
strings are topologically stable objects which might be found during a phase 
transition in the early universe (Kibble [1]). Cosmic strings play an important 
role in the study of the early universe. These arise during the phase transition 
after the big bang explosion as the temperature goes down below some critical 
temperature as predicted by grand unified theories (Zel'dovich et al. [2]; Kibble 
[1, 3]; Everett [4]; Vilenkin [5]). It is believed that cosmic strings give rise 
to density perturbations which lead to the formation of galaxies (Zel'dovich 
[6]). These cosmic strings have stress-energy and couple to the gravitational 
field. Therefore it is interesting to study the gravitational effects that arise 
from strings. 
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The general relativistic treatment of strings was initiated by Letelier [7, 8] 
and Stachel [9] . Letelier [7] has obtained the solution to Einstein's field equations 
for a cloud of strings with spherical, plane and cylindrical symmetry. Then, in 
1983, he solved Einstein's field equations for a cloud of massive strings and 
obtained cosmological models in Bianchi I and Kantowski-Sachs space-times. 
Benerjee et al. [10] have investigated an axially symmetric Bianchi type I string 
dust cosmological model in presence and absence of magnetic field. The string 
cosmological models with a magnetic field are also discussed by Chakraborty 
[11], Tikekar and Patel [12]. Patel and Maharaj [13] investigated stationary 
rotating world model with magnetic field. Ram and Singh [14] obtained some 
new exact solutions of string cosmology with and without a source free magnetic 
field for Bianchi type I space-time in the different basic form considered by 
Carminati and Mcintosh [15]. Exact solutions of string cosmology for Bianchi 
type II, VIq, VIII and IX space-times have been studied by Krori et al. [16] and 
Wang [17]. Singh and Singh [18] investigated string cosmological models with 
magnetic field in the context of space-time with G3 symmetry. Lidsey, Wands 
and Copeland [19] have reviewed aspects of super string cosmology with the 
emphasis on the cosmological implications of duality symmetries in the theory. 
Baysal et al. [20] have investigated the behaviour of a string in the cylindrically 
symmetric inhomogeneous universe. Bali et al. [21 - 25] have obtained Bianchi 
types IX, type-V and type-I string cosmological models in general relativity. 
Yavuz [26] have examined charged strange quark matter attached to the string 
cloud in the spherical symmetric space-time admitting one-parameter group of 
conformal motion. Recently Kaluza-Klein cosmological solutions are obtained 
by Yilmaz [27] for quark matter coupled to the string cloud in the context of 
general relativity. 

On the other hand, the matter distribution is satisfactorily described by per- 
fect fluids due to the large scale distribution of galaxies in our universe. How- 
ever, a realistic treatment of the problem requires the consideration of material 
distribution other than the perfect fluid. It is well known that when neutrino 
decoupling occurs, the matter behaves as a viscous fluid in an early stage of the 
universe. Viscous fluid cosmological models of early universe have been widely 
discussed in the literature. 

Recently, Seljak et al. [28] have studied cosmological parameters from com- 
bining Lyman-a forest with C.M.B. They established the cosmic strings limited 
to Gfi < 2.3 x 10~ 7 at 95 percent cathodoluminescence (c.l.), and correlated iso- 
curvature models are also tightly constrained. Observations on C.M.B. constrain 
the cosmic string which is predicted in terms of cosmic string power spectrum 
[29]. Howerever, we are not concerned with this new phenomena. Tikekar and 
Patel [12], following the techniques used by Letelier and Stachel, obtained some 
exact Bianchi III cosmological solutions of massive strings in presence of mag- 
netic field. Maharaj et al. [30] have generalized the previous solutions obtained 
by Tikekar and Patel [12] considering Lie point symmetries. Yadav et al. [31] 
have studied some Bianchi type I viscous fluid string cosmological models with 
magnetic field. Recently Wang [32 - 35] has also discussed LRS Bianchi type I 
and Bianchi type III cosmological models for a cloud string with bulk viscos- 
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ity. Very recently, Yadav, Rai and Pradhan [36] have found the integrability of 
cosmic string in Bianchi type III space-time in presence of bulk viscous fluid by 
applying a new technique. Motivated by the situations discussed above, in this 
study, we established a formalism for studying the new integrability of massive 
strings in Bianchi III space-time in presence of a variable bulk viscosity. 



2 The Metric and Field Equations 

We consider the space-time of general Bianchi III type with the metric 

ds 2 = -dt 2 + A 2 [t)dx 2 + B 2 {t)e~ 2ax dy 2 + C 2 (t)dz 2 , (1) 

where a is constant. The energy momentum tensor for a cloud of string dust 
with a bulk viscous fluid of string is given by Letelier [7] and Landau & Lifshitz 
[37] 

T'l = pr,r ! - \ Xi x j - fuj, (g> + r,r ! ) , (2) 

where Vi and X{ satisfy condition 

v % vi = -x l x l = — 1, v l Xi = 0, (3) 

where p is the proper energy density for a cloud string with particles attached to 
them, A is the string tension density, v l is the four-velocity of the particles, and 
x % is a unit space-like vector representing the direction of string. If the particle 
density of the configuration is denoted by p p , then we have 

p = Pp + X. (4) 

The Einstein's field equations (in gravitational units c — 1, G — 1) read 

R{ -\R9l= -8*2* (5) 

where Rj is the Ricci tensor; R = g v Rij is the Ricci scalar. In a co-moving 
co-ordinate system, we have 

^ = (0,0,0,1), x* = (0,0,1/C,0). (6) 

The field equations ([5]) with Eq. ||5J) subsequently lead to the following system 
of equations: 

B 44 C 44 B 4 C 4 

-lT + — + W =8 ^' (?) 

A44 A4C4 

~A + ~C + AC 

A 44 -B44 A 4 B 4 a 2 
~T + ~B~ + AB ~ A 2 " 

AaBa BaCa CaAa a 2 , . 

-AB- + ^C + ^A-A^= 8 ^ (10) 



8tt#, (8) 
= 87r(A + ^), (9) 
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0, 



(11) 



Aa_Ba 
A B 

where the subscript 4 at the symbols A, B and C denotes ordinary differentiation 
with respect to t. The particle density p p is given by 



8np p 
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(12) 



in accordance with Eq. (j4]). 



The velocity field v % as specified by Eq. ^ is irrotational. The expansion 
(9) and components of shear tensor (<7y) are given by 
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(13) 

(14) 

(15) 

(16) 
(17) 



(18) 



3 Solutions of the Field Equations 

The field equations (|7)l- (fTT|) are a system of five equations with six unknown 
parameters A, B, C, p, A and £. One additional constraint relating these param- 
eters is required to obtain explicit solutions of the system. Referring to Thorn 
[38] , observations of the velocity-redshift relation for extragalactic sources sug- 
gest that Hubble expansion of the universe is isotropic today to within about 
30 percent [39, 40]. Put more precisely, redshift studies place the limit 

I < 0-30 

on the ratio of shear, a, to Hubble constant H in the neighbourhood of our 
Galaxy. Following Bali and Jain [41] and Pradhan et al. [42, 43], we assume 



4 



that the expansion (9) in the model is proportional to the shear (tr), which is 
physically plausible condition. This condition leads to 

B = £C n (19) 

where I is proportionality constant and n is a constant. Equations leads to 

A = mB, (20) 

where m is an integrating constant. To obtain the determinate model of the 
universe, we assume that the coefficient of bulk viscosity is inversely proportional 
to expansion (6). This condition leads to 

= k, (21) 

where K is a proportionality constant. 

By the use of Eqs. (JT^ and $ZT$, Eq. (7J) leads to 

2n 2 \ C\ (IQitK\ 



Let 

Ci = /(C). (23) 

Hence Eq. ([22]) leads to 

From Eq. (pM)) , we have 

f*C& = (^»^)c !s£ &* + L i (25) 



y .i 2 + n+l / 

where L is an integrating constant. Eq. (|25[) leads to 



/ 2 =(^^T)C 2 + LC-W. (26) 
1 n z + n + 1 ' 



From Eq. (j26|) , we have 

C^dC 



sfc 2 J^^) +b 2 
Integrating Eq. ([27]) . we have 



y/adt. (27) 



C 



f3 2 (pt + 1 ) 2 -b 2 
1 



(28) 
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where 7 is an integrating constant and 

{n 2 +n + l)y/a 



P 



Therefore, we have 

B = ec n = £ 

A = mB = £m 



2(n+l) 

8ttK 
(n 2 + n + l))' 

.0 L 



p 2 {pt + 1 ) 2 -b 2 



P 2 {pt + 1 ) 2 -b 2 
P 2 



( n(n+l) \ 



( ~(n+l) \ 



Hence the metric ([1} reduces to the form 



ds 2 = -dt 2 



P 2 {pt + 1 ) 2 -b 2 
P 2 



[dX 2 + , 



- 2aX dY 2 ] 



[3 2 {f3t + 1 ) 2 -b 
P 2 



2(" + l) 



dZ 2 



where 



imx = X, £y — Y, z = Z. 

Using the transformation 

Pifit+i) = 6cos(/3r) 
the metric (|3"Tj) reduces to the form 

2 _ & 2 sin 2 Q3t) 2 [ 6 4 sin 4 (/3r) ] ^TT 
/3 2 + [ /? 4 

6 4 sin 4 (/3r)l^TT^ 2 



(29) 
(30) 



(31) 



(32) 



dX 2 + e—^dY 2 



P 4 



(33) 



In absence of viscosity i.e. when j3 — > then the metric (f3"3")l reduces to the new 
form 



ds 2 = -b 2 r 2 dr 2 + (6 4 r 4 ^ — 



(6 4 r 4 )^^rfZ 2 . (34) 
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4 Some Physical and Geometrical Features of 
the Models 

The rest energy density (p), the string tension density (A), the particle density 
(p p ), the scalar of expansion (#) and the shear (a) for the model (|31 [) are given 

by 



8irp 



An{n + 2)(n+l) 2 /3 2 (/?t + 7 ) 



8ttA 
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(36) 
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/3 2 



(n 2 + n + 1) 
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/9 2 



/3 2 (n + l) 2 
(n 2 +n + 1)' 



2(2n + l)(n + l)/3(/ft + 7 ) 



2(n 2 - l)/3(/3t + 7) 



(n 2 + n + 1)V3 



/3 2 (/3t+7) 2_ b 2 



(37) 
(38) 
(39) 



The energy conditions p > and p p > are satisfied in the presence of bulk 
viscosity for the model (|3"Tj) . The condition p > leads to 



a(n 2 +71+1) 



2£m(n + i)y/n(n + 2) 



< 



< 



a(n 2 + n + 1) 



3" 



2fm(n + l)^/n(n + 2) 



From Eq. (j36|) . we observe that the string tension density A > provided 

1 [fit + 7) 2 (n 2 +71 - 2) 
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The model in presence of bulk viscosity starts with a big bang at time t = -Jr — 5 . 
The expansion in the model decreases as time increases. The expansion in the 
model stops at t = — The model ([3"Tj) in general represents shearing and non- 
rotating universe. The role of bulk viscosity is to retard expansion in the model. 
We can see from the above discussion that the bulk viscosity plays a significant 
role in the evolution of the universe. Furthermore, since limt^oo § 7^ 0, the 
model does not approach isotropy for large value of t. However the model 
isotropizes when t = — ^ . There is a real physical singularity in the model (|31[) 
at t = 0. 

Using the transformation 

+ 7) = 6cos(/3r), 

the rest energy density (p), the string tension density (A), the particle density 
(p p ), the scalar of expansion (9) and the shear (er) for the model (|33[) are given 

by 



8np- 



An{n + 2)(n + I) 2 cos 2 (/3r) 



b 2 (j 
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8np p 
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where b = —N, N > 0. In absence of bulk viscosity i.e. when (3 - 
physical and kinematic quantities for the model (|34[) are given by 
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S^Ll (48) 



N(n 2 + n + 1) t 

2(n 2 - 1) 1 
7V\/3(n 2 + n + 1) r 
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(49) 



In the absence of bulk viscosity, the energy conditions p > and p p > are 
satisfied for the model (|3~4")) . The condition p > leads to 

4n(n + l) 2 (n + 2)to 2 £ 2 0™ 2 +™- 1 ) 



rn 2+ n +l < 



a 2 (n 2 + n + l) 2 
The string tension density A > if 

An(n + l)(n 2 + n - 2)m 2 £ 2 b 2( - n2+n -^ 



T „2 + n + 1 < 



a 2 (n 2 + n + l) s 



The model in the absence of bulk viscosity starts with a big bang at 
time r = and the expansion in the model decreases as time increases. Since 
limT-^oo 5 =^ 0, the model does not approach isotropy for large value of r. 
However the model isotropizes for n = — 1 and n = 1. There is a real physical 
singularity in the model ([34]) at r = 0. 
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